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Abstract 

Diurnal  variation  is  one  of  the  significant  instabilities  in  Two-Way  Satellite  Time  and 
Frequency  Transfer  (TWSTFT)  data.  It  is  necessary  to  correct  the  diurnal  variation  for 
comparing  the  time-scale  difference.  We  focus  on  the  up-ldownlink  delay  difference  caused  by 
satellite  motion.  In  this  paper,  we  propose  to  correct  the  TWSTFT  data  by  using  round-trip 
delay  measurement.  There  are  four  laboratories:  the  National  Institute  of  Information  and 
Communications  (NICT),  Korea  Research  Institute  of  Standards  and  Science  (KRISS),  the 
National  Time  Service  Center  (NTSC),  and  TL,  who  join  the  TWSTFT  network  by  using  the 
NICT  modem.  Moreover,  TL  and  NICT  are  continuing  the  DPN  experiment.  In  the  end,  we 
investigate  these  TWSTFT  data  and  the  corrected  data  are  performed. 


INTRODUCTION 

TWSTFT  is  one  of  the  most  precise  technologies  for  comparing  time-scale  differences.  Because 
two-way  signal  paths  are  highly  reciprocal,  it  can  achieve  precise  stability.  There  is  another  method, 
GPS  Precise  Point  Positioning  (GPS  PPP),  for  comparing  time-scale  differences.  However,  GPS  PPP’s 
data  need  to  be  arranged  by  post-processing.  It  is  an  advantage  that  TWSTFT  can  perform  time-scale 
differencing  in  real  time. 

However,  TWSTFT  faces  a  serious  problem.  Diurnal  variation  is  one  of  the  sources  of  instability  and 
we  find  it  significant  in  TWSTFT  results.  Since  we  do  not  find  the  diurnal  from  GPS  PPP’s  results,  it  is 
impossible  to  evaluate  the  stability  of  time-scale  differences  by  using  TWSTFT.  Many  researchers  work 
at  finding  the  diurnal  sources.  [1]  mentioned  the  variation  of  Sagnac  residual  according  to  ephemeris 
data,  and  corrected  TWSTFT  results  for  it.  [2]  proposed  a  correction  method  by  adopting  the 
troposphere  delay  model.  [3]  indicated  and  analyzed  the  possible  instabilities,  including  environmental 
factors,  ionosphere  delay,  satellite  motion,  and  the  satellite  transponder.  [4]  further  investigated  the 
influence  of  ionosphere  effect,  satellite  motion,  and  environmental  conditions. 

It  is  well-known  that  in  order  to  keep  the  satellite  static  along  the  geostationary  orbit,  the  satellite’s 
operator  has  to  adjust  the  position  of  satellite.  The  satellite  will  be  set  in  slow  motion  and,  thus,  have  a 
relative  velocity  with  respect  to  the  earth  station.  It  causes  the  up-/downlink  signal  paths  to  be  no  longer 
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the  same.  That  is  called  the  integrated  Doppler  shift  [5].  This  paper  is  concerned  with  the  correction 
based  on  the  integrated  Doppler  shift.  In  this  paper,  we  propose  a  method  to  calculate  the  integrated 
Doppler  shift  and  then  correct  the  TWSTFT  results. 


TWSTFT  NON-RECIPROCITY 

The  following  acronyms  are  defined  according  to  ITU-R  TF.  1153-2  [6].  In  this  paper,  we  denote  the 
time-scale  difference  between  I PPSTXx  and  \PPSTX2  as  dt ,  and  dt  is  either  greater  or  less  than 

zero.  By  the  TWSTFT  principle  depicted  in  Fig.  1,  at  station  1,  the  1PPS  signal  from  station  2  will  be 
received  and  the  following  measurement  will  be  shown  on  the  time-interval  counter. 

77,  =dt  +  TX2  +  SPU2  +  SCU2  +  SPT2  +  SPDX  +  SCO,  +RXt  (1) 

The  1PPS  signal  from  station  1  will  be  measured  and  shown  at  station  2. 

77,  =  -dt  +  TXx+  SPUt  +  SCUX  +  SPT  +  SPD2  +  SCD2  +  RX,  (2) 

For  simplicity,  we  combine  the  up-/downlink  delay  with  the  Sagnac  correction.  Let 
SDl=SPDx+SCDl  -  then 


77,  =  dt  +  TX2  +  SU2  +  SPT2  +  SDl  +  RXx  (3) 

and 

TI2  =  -dt  +  TXx+  SUX  +  SPTX  +  SD2  +  RX2  (4). 

The  time-scale  difference  is  given  by  adding  (3)  with  (4). 

dt  =  0.5  x  [77,-77,] 

+  0.5  x  {[517, -5D,]-[5I/2-5D2]} 

+  0.5  x[SPTx-SPT2] 

In  the  right-hand  side  of  (5),  the  first  term,  0.5  x  [77,-77,],  represents  the  difference  of  counter 

readings,  which  we  call  TWSTFT  data  in  this  paper.  The  other  terms  address  non-reciprocal  factors. 
In  this  paper,  we  focus  on  the  up-/downlink  difference  corresponding  to  the  third  term. 
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Fig.  1.  Conceptual  representation  of  time-interval  counter  readings,  77,  and  772 , 
along  the  two-way  signal  path. 


Fig.  2.  Conceptual  representation  of  time-interval  counter  readings  (77,R  and  77,s ) 
along  the  round-trip  signal  path. 


NOTATION  AND  ASSUMPTION 

Denote  v,  as  the  rate  of  change  of  the  distance  R;  (/)  between  earth  station  i  and  the  satellite.  In 
other  words,  v,  is  regarded  as  the  relative  velocity,  and  has  the  following  relationship  with  the  relative 
distance  /?(/). 

Ri(t)  =  Ri(0)-vi-t  (6) 

Here,  we  assume  vi  is  fixed  during  300  seconds,  and  has  a  plus  sign  when  the  satellite  approaches  the 
earth  station.  Because  station  1  has  to  receive  its  1PPS  signal  along  the  round-trip  path,  we  further 
denote  the  round-trip  downlink  delay  as  SD ,R ,  different  from  SDl ,  which  is  the  downlink  delay  along 
the  two-way  path.  Figure  2  depicts  the  delay  of  each  segment  along  the  round-trip  path.  We  ignore  the 
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reference  delay  among  all  laboratories  and  the  ionosphere  delay  at  each  site,  because  those  have  diurnal 
variations  within  hundreds  of  picoseconds,  while  the  up-/downlink  delay  would  have  them  in  the 
hundreds  of  microseconds.  We  also  ignore  the  transponder  delay  because  that  is  not  known.  We 
further  assume  the  1PPS  signal  travels  in  free  space;  thus,  the  signal  speed,  defined  as  c ,  is  a  constant 
equal  to  299,792,458  meters/second. 


INTEGRATED  DOPPLER  SHIFT 

This  section  is  aimed  at  deriving  the  up-/downlink  non-reciprocity  caused  by  satellite  motion.  If  the 
satellite  is  strictly  geostationary,  then  the  uplink  and  downlink  delays  are  the  same  when  the  1PPS  signal 
travels  along  the  round-trip  path.  That  is, 


p  R,(  0) 

SU,  =  SD f  =  (7). 

c 

What  happens  if  the  satellite  is  in  motion?  First,  we  consider  the  1PPS  signal  transmitted  from  station  1. 
We  can  get  uplink  delay  SUl  from  (6). 


SUl  c  =  R{  (0)-5C/,  ■vl  (8) 

The  satellite  will  receive  and  bypass  the  1PPS  signal  from  station  1  after  SUl  seconds.  Along  the 
round-trip  path,  the  downlink  delay  is,  thus,  given  by  (6). 

SD*  ■c  =  Ri  (5C7, )  -  SD*  ■  v,  (9) 

Along  the  two-way  path,  the  downlink  delay  is  also  given  by  (6). 

SD2  ■c  =  R2(SUt)- SD2  ■  v2  (10) 

From  (8)  to  (10),  we  can  get  the  propagation  delay  of  each  segment  in  terms  of  relative  distance  and 
velocity. 


SU,  = - R,(  0)  (11) 


c  +  v, 


SDIg=^-^TRl(  0)  (12) 


(c  +  v,)- 


SD2  = 


c  +  v. 


(c  +  V,)(c  +  V2) 


Rt  (o)  (13) 


On  the  other  hand,  we  can  get  the  propagation  delay  when  we  consider  the  1PPS  signal  transmitted  from 
station  2. 
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SU2  = 


c  +  v. 


( C  +  V2 ) 


1 


SD,= - R,(d t)— - -j - - 

c  +  Vi  (c  +  v1)(c  +  v2) 


R2(dt)  (14) 
R2(dt)  (15) 
R2(dt) 


(16) 


It  is  obvious  that  the  propagation  delays  in  (11)  and  (12)  are  different  from  those  in  (7).  The  difference 
caused  by  satellite  motion  is  the  so-called  integrated  Doppler  shift.  Therefore,  from  (11)  to  (16),  the 
up-/downlink  non-reciprocity  that  appears  in  the  third  term  of  (5)  can  be  expressed  in  terms  of  relative 
distance  and  velocity. 


0.5  x  {[SUl  -  SD,  ]  -  [SU2  -  SD2  ]} 


V2^1  (0)— V1  R2{dt) 

2(c  +  Vi)(c  +  v2) 


(17) 


(17)  is  called  the  integrated  Doppler  correction.  To  derive  (17),  we  assume  Rx  (dt )  =  Rt  (0)  —  dt  ■  v1  is 
approximately  equal  to  R^  (0) .  It  is  reasonable  that  the  satellite  does  not  move  rapidly  within  a  few 
tenths  of  nanoseconds.  In  principle,  dt  •  Vj  is  at  the  level  of  a  few  tenths  of  nanometers,  and  can  be 
ignored  if  R{  (0)  is  about  40,000  kilometers.  In  the  next  section,  we  will  use  the  round-trip  delay 
measurement  to  estimate  the  relative  distances  Rt  ( 0 )  and  R2  (dt  j  ,  and  relative  velocities  Vj  and  v2. 


ESTIMATION  USING  ROUND-TRIP  MEASUREMENT 


The  round-trip  signal  path  is  the  path  along  which  the  1PPS  signal  travels  and  returns  between  the  earth 
station  and  the  satellite.  Only  station  1  is  considered  in  this  section.  From  (8)  and  (9),  the  round-trip 
delay  can  be  expressed  as  follows. 


77,*(0)  :=SUt  +  SB*  =  (2C  +  V'J  Rt(Q)  (18) 

(c+vq) 


We  define  TId  (k  ')  as  the  round-trip  delay  measurement  at  time  k  when  beginning  to  transmit  the 
1PPS  signal  at  station  1.  From  (6),  the  round-trip  delay  measurement  yields 


(c  +  Vi)  (c  +  Vj)  (c  +  Vj) 


k  =  a-bk  (19). 


Assume  that  (19)  is  a  signal  model  and  we  observe  TI*  (k)  for  k  =0,1,..., N  —  1.  Then,  according  to 
the  least-squares  approach,  we  can  estimate  a  and  b  by  minimizing  the  following  equation. 
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/(«,/?)  =  jr(77  ?(k)-(a-b-kjf  (20) 


k= 0 


By  differentiating  with  respect  to  a  and  b  and  setting  the  results  to  zero,  we  get 

f  N-l 


a  = 


b  = 


B-Y,TI?(k)-A-Y,(lc-TI?(k))\  (B-A2) 

V  k= 0  k= 0  )t 


f  N-l 


N-l 


where  < 


A.^Tlf^-^k-TI^k))  ( B-A 2) 


V  k= 0 


k= 0 


A  = 


B  = 


\n-l) 


(21). 


~  2cv,  +  v,2 

From  (19)  and  (21),  we  set  b  = - —  and  then  we  can  get  the  estimated  relative  velocity. 

0  +  v.) 


l-i/l-i 

V,  = — ,  c  (22) 


'  T 


-b 


2  c  +  v, 


In  the  same  way,  we  set  d  = - x—  R]  ( 0 )  and  then  we  can  get  the  estimated  relative  distance. 

(c  +  v  i) 

,  ,  (c  +  V,)2  „ 

(o)  =  - —a  (23) 

2c  + Vj 


We  can  follow  the  same  procedure  to  derive  v2  and  R2  ( dt )  . 


DOPPLER  EFFECT 

In  this  section,  we  will  verify  that  our  derivation  is  identical  to  the  Doppler  Effect.  If  the  satellite  has  a 
relative  velocity  with  respect  to  an  earth  station,  we  will  find  the  signal  frequency  observed  at  the  satellite 
changed.  The  following  formula  is  the  Doppler  Effect  regarding  all  kinds  of  electromagnetic  waves 
traveling  in  the  free  space: 


4f  =  v 
/  c 


(24) 


where  we  denote  /  ,  A/  and  v  as  the  signal  frequency,  the  Doppler  frequency,  and  the  relative 
velocity.  From  this  formula,  we  can  get  the  frequency  fSAT  observed  at  the  satellite  when  station  1 
transmits  a  1PPS  signal  with  frequency  fx . 

(c  +  v, ) 

Asat  =  - - -fi  (25) 

c 


The  1PPS  signal  will  be  bypassed  and  reflected  back  to  station  1  by  the  satellite,  and  we  can  measure  the 
frequency  round_trip  along  the  round-trip  path  at  station  1 . 
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(c  +  V.)  (c  +  V,)2 

=  ' - -  Asat  =  /i  (26) 


From  this  point  of  view,  we  verify  our  derivation  by  giving  a  short  example.  Since  1  second  is  the  time 
interval  of  two  consecutive  1PPS  signals,  by  transmitting  two  consecutive  1PPS  signals,  we  can  receive 
them  along  the  round-trip  path  and  measure  the  time  interval  T'  between  them.  From  (19),  we  have 


1,  round  —trip 


.  =i  +  27*(l)-r/*(0)  = 


(2 c  +  v,)  ,  ,  (2c  +  v, ) 

1+2 - 4rx{  i)-2 - 4 


(c+vi  y 


(c  +  Vi)‘ 


(0)  = 


(c  +  Vi)- 


(27). 


Equivalently,  the  observed  frequency  of  1PPS  signal  will  be  the  inverse  of  the  time  interval  T ' .  That  is, 

1  _(c  +  v1)2 


A 


,  round -trip 


(28). 


1 1,  round -trip 


This  is  the  same  as  (26)  when  f\  is  set  to  1  Hz,  the  frequency  of  1PPS  signal. 


EXPERIMENT  RESULTS 

We  have  to  receive  the  signals  along  the  two-way  and  round-trip  path  at  the  same  epoch,  so  it  is  necessary 
to  use  the  multi-channel  modem  to  get  them.  At  present,  the  NICT  modem  is  a  multi-channel  modem 
and  is  being  employed  among  the  Asia-Pacific  TWSTFT  links.  Besides,  Gotoh  and  Amagai  [7]  have 
developed  a  novel  technique  of  TWSTFT  with  high  precision.  They  use  Dual  PRN  codes  (DPN)  signal 
as  the  1PPS  signal  and  design  a  GPU-based  receiver  by  parallel  signal  processing.  TL  and  NICT  are 
now  continuing  the  DPN  experiment  by  which  a  round-trip  delay  measurement  is  also  available.  Both 
experiments  are  performing  via  the  geostationary  satellite  IS -8  at  166E,  and  those  data  exchange  though 
the  Internet. 

We  regard  the  first  terms  in  the  right-hand  side  of  (5)  as  TWSTFT  data.  The  TWSTFT  data  from  NICT 
modem  and  from  DPN  are  depicted  in  Fig.  3.  In  Fig.  3,  we  can  find  that  the  diurnal  variation  exists  in 
the  NICT  and  in  the  DPN  modem,  and  their  tendencies  are  the  same.  However,  their  magnitudes  of 
diurnal  variation  are  not  the  same.  It  is  evident  that  we  have  to  use  different  multiple  factors  for 
correcting  TWSTFT  data.  Practically,  a  factor  of  (17)  will  be  used  for  correcting  TWSTFT  data. 

The  TWSTFT  data  and  time  deviation  (TDEV)  results  are  depicted  in  Figs.  4  to  6.  The  2-month  results 
are  shown.  We  use  (17)  by  a  factor  of  40  to  correct  TWSTFT  data.  Obviously,  diurnal  variation 
decreases,  but  does  not  vanish  after  TWSTFT  data  are  corrected.  There  would  be  another  source  of 
instability  that  causes  diurnal  variation. 
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Fig.  3.  TWSTFT  data  of  NICT-TL  link  from  the  NICT  modem  (red  line)  and  the  DPN 
modem  (green  line)  with  averaging  times  of  1  second  (left)  and  300  seconds  (right). 
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Fig.  4.  TWSTFT  data  (left)  and  TDEV  (right)  of  the  NICT-KRISS  link  from  the  NICT  modem. 
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NICT  -  NTSC  NICT  -  NTSC ,  TDEV 


Fig.  5.  TWSTFT  data  (left)  and  TDEV  (right)  of  the  NICT-NTSC  link  from  the  NICT  modem. 


NICT-TL  NICT  -  TL.  TDEV 


Fig.  6.  TWSTFT  data  (left)  and  TDEV  (right)  of  the  NICT-TL  link  from  the  NICT  modem. 
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NICT  -  TL  (DPN)  NICT  -  TL  (DPN) ,  TDEV 


Fig.  7.  TWSTFT  data  (left)  and  TDEV  (right)  of  the  NICT-TL  link  from  the  DPN  modem. 


Also,  we  show  the  7-day  TWSTFT  data  from  the  DPN  modem  in  Fig  7.  We  find  that  the  diurnal 
variation  decreases  with  TWSTFT  data  corrected  by  a  factor  of  20. 


Round-Trip  Doppler  Frequency  @  TL 


MJD 

Fig.  8.  Measured  Doppler  frequency  (red  line)  and  Doppler  frequency  derived  by  (24)-(26)  (green  line). 
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We  verify  that  our  derivation  of  relative  velocity  matches  the  theoretical  Doppler  frequency.  Note  that 
the  transmit  frequency  is  14.164  GHz,  the  change  of  frequency  is  about  1753  MHz  due  to  of  transponder, 
and  the  receive  frequency  is  12.41 1  GHz.  The  measured  frequency  is  recorded  by  the  DPN  modem.  In 
Fig.  8,  we  can  find  that  the  two  curves  are  in  good  agreement;  thus,  our  derivation  is  meaningful  and 
practical. 


CONCLUSION 

We  use  the  integrated  Doppler  correction  to  improve  TWSTFT  stability,  and  the  diurnal  variation  drops 
considerably.  We  consider  that  the  integrated  Doppler  shift  may  be  the  dominant  factor  of  diurnal 
variation  even  if  diurnal  variation  does  not  vanish  after  correction.  From  experimental  results,  we  find 
that  the  theoretical  Doppler  frequency  is  in  good  agreement  with  the  measured  Doppler  frequency.  We 
also  find  that  integrated  Doppler  correction  needs  a  multiple  factor  because  the  magnitudes  of  diurnal 
variation  are  not  the  same  in  the  NICT  and  the  DPN  modems.  It  may  be  caused  by  signal  waveform 
distortion,  which  we  will  investigate  in  the  near  future. 
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